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Abstract Understanding the inherent damping mechanisms of floor vibrations has 
become a matter of increasing importance following the development of new 
composite floor layouts and increased span. The present study focuses on the 
evaluation of material damping in timber beam specimens with dimensions that are 
typical of common timber floor structures. Using the impact test method, 11 solid 
wood beams and 11 glulam beams made out of Norway Spruce (Picea abies) were 
subjected to flexural vibrations. The tests involved different spans and orientations. 
A total of 420 material damping evaluations were performed, and the results are 
presented as mean values for each configuration along with important statistical 
indicators to quantify their reliability. The consistency of the experimental method 
was validated with respect to repeatability and reproducibility. General trends found 
an increasing damping ratio for higher modes, shorter spans, and edgewise orien¬ 
tations. It is concluded from the results that material damping of timber beams of 
structural dimensions is governed by shear deformation, which can be expressed 
more conveniently with respect to the specific mode shape and its derivatives. 


Introduction 


Damping is of interest in many disciplines and was first investigated with regard to 
space structure applications (Neumark 1962). Nowadays, it is also relevant for more 
common structures, for example, cars or buildings, under the generic name of 
“comfort properties.” When designing building floors, excessive vibrations are 
usually not a safety concern, rather a serviceability issue due to annoyance and other 
discomforts. When damping is close to resonance, it will have a large beneficial 
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influence on the structural response since it decreases both the amplitude of steady- 
state oscillations as well as the duration of transient oscillations. Despite its 
substantial effects, damping is rarely prescribed in design codes or standards, 
mainly because of the lack of knowledge. For timber structures, damping 
evaluations are generally considered to depend too much on the engineer’s 
judgment, because of the lack of reliability in the methods. Moreover, the poor 
understanding of what the damping phenomenon is exactly leads to the situation in 
Eurocode 5 (European Committee for Standardization 2005): Damping quantity is 
introduced as an arbitrary single value, nationally determined for any type of floor. 
A cruder criterion is commonly used that employs the natural frequency as a target 
value to assess the performance of the structure with respect to its comfort 
properties. This is despite the fact that the implementation of damping in prediction 
models would result in more accurate predictions and more cost-effective building 
solutions. 

Data on the measured damping properties of wood are scarce, partly due to 
difficulties in retrieving reliable experimental data. In addition, damping nomen¬ 
clature is often inconsistent (Lazan 1968), which increases the confusion by mixing 
different types and quantities of damping. Among others, the damping ratio £, the 
logarithmic decrement 2, the internal friction g -1 , the loss factor rj , or the loss 
tangent tan(^) were successively used in different studies. They are all related by: 

A Q~ l r\ tan(<3) 

271 2 2 2 

From a material science point of view, material damping may be defined as 
internal friction by transformation of mechanical energy into heat during cyclic 
stress (Ouis 2002). In general, damping may be divided into three classes 
(Woodhouse 1998): 

• internal friction throughout the material making up the structure: material 

damping, 

• energy dissipation associated with junctions or interfaces between parts of the 

structure: structural damping, and 

• energy dissipation associated with a fluid in contact with the structure, involving 

either local viscous effects or radiation into the fluid: fluid damping. 

The making of musical instruments seems to have induced scientific investiga¬ 
tions of wood properties, including damping estimations, much earlier than 
corresponding investigations of timber structures. In order to improve knowledge on 
wood for musical instruments, the frequency dependence of the logarithmic 
decrement was investigated (Fukada 1950), as well as the effects of moisture 
content and temperature (Fukada 1951). Particular attention was paid to the 
selection of suitable types of wood for the soundboards of musical instruments, and 
relationships between internal friction and dynamic Young’s modulus were 
investigated to that end (Ono and Norimoto 1985). Investigations about the 
dependence of internal friction to frequency were numerous for wood specimens 
used in musical instruments, but were not always concordant. On the one hand, 
damping was found to be independent of frequency within large ranges of 
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frequency, for example, up to about 2 kHz (Ono and Norimoto 1985) or up to 
10 kHz (Nakao et al. 1985). On the other hand, damping was observed to be 
independent of frequency (Foster 1992), but only in the longitudinal direction, 
whereas it was strongly dependent on frequency in the radial direction (Foster 
1992). The influence of the grain angle on the vibration properties, including 
damping, was also investigated more generally (Bremaud et al. 2011). 

Experimental damping values have been largely available for an acoustic use of 
timber (Obataya et al. 2000; Ono and Norimoto 1985; Bucur 2006; Spycher et al. 
2008), or else concerning the mechanical pulping use of timber (Havimo 2009). 
More recently, 54 tropical species were analyzed (Brancheriau et al. 2010) using 
free vibration and forced-released vibration tests on 350-mm-long specimens. 

However, few studies in the literature report on the damping properties of timber 
elements for structural use. Even if, as early as 1927, Kimball and Lovell (1927) used 
one-meter-long specimens to measure the internal friction of 18 different solids 
including maple wood, and experimental damping values related to a structural use of 
timber were first reported by Yeh et al. (1971). They performed free flexural vibration 
tests on real-sized wood frames to investigate both the damping associated with the 
material itself and the damping associated with the joining device, such as nails or glue. 

Due to the scarce quantification of damping properties for structural use, the 
intention of the present study is to provide new and reliable values for material 
damping of timber beams that are typical for common floor structures. This work 
first describes an experimental method to evaluate material damping in timber 
beams of structural dimensions. The statistical methods used for evaluating 
reliability are defined in “Materials and methods.” In “Result and discussion,” the 
validity of the method is assessed with respect to repeatability and reproducibility, 
and the reliability of the results is demonstrated by means of statistical indicators. 
Governing parameters are also discussed. Some of the cited experimental damping 
measurements are reported in Table 1 for ease of comparison. Damping is expressed 
as damping ratio values £, as defined in Eq. (1). 


Materials and methods 


Experimental setup 


A total of 22 beams were tested using flexural vibration estimation techniques. 
These were 11 solid wood beams of quality C24 (European Committee for 
Standardization 2009) and 11 glulam beams of quality GL32c (European 
Committee for Standardization 1999), all made out of Norway Spruce (Picea 
abies). The nominal dimensions and properties are summarized in Table 2, and they 
represent typical dimensions in common timber floor structures. All beams were 
measured and weighed to determine the density. Glulam beams had a mean density 

Q 

of 494 kg/m (coefficient of variation = 2 %), while the solid wood beams had a 

Q 

mean density of 445 kg/m (coefficient of variation = 6 %), both match the 
nominal properties of their respective strength class. The orientation of annual rings 
was not recorded. 
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Table 1 Experimental damping ratio presented in previous studies 


Author(s) 


Damping ratio Wood species 


Method 




Fukada (1950) 

Kollmann and 

Krech (1960) 

Matsumoto 

(1962) 

Holtz cited in 
Bucur (2006) 

Yeh et al. (1971) 

Wert et al. (1984) 


0.0027 

0.0031 


Yezo Spruce ( Picea Jezoensis) 
Norway Spruce {Picea abies) 


Free flexural vibration 
Forced flexural vibration 


0.0027 


Japanese Cedar ( Cryptomeria 
japonica ) 

Spruce {Picea abies ) 


Transverse vibration 


0.0030 


0.0035 

0.0075 


Hemlock ( Tsuga ) 

Sitka Spruce {Picea sitchensis) 


Free flexural vibration 

Fow inverted torsional 
pendulum 

Free longitudinal vibration 


0.0056 

(longitudinal) 


Ono and 
Norimoto 

(1985) 

Nakao et al. 

(1985) 

Foster (1992) 


Sitka spruce {Picea sitchensis ) 


0.0035 


Sitka spruce {Picea sitchensis) 


Free flexural vibration 


0.0030 

(axially) 

0.0025 

(axially) 

0.0022-0.0037 


King William Pine 

{Athrotaxis selaginoides ) Norway 
Spruce {Picea abies ) 


Uniaxial loading + phase 
lag measurement 


Obataya et al. 

2000) 

Spycher et al. 
(2008) 


Sitka Spruce {Picea sitchensis ) 


Free flexural vibration 


0.0051 

(axially) 

0.0068 

(axially) 

0.0035 


Norway Spruce {Picea abies ) 
Sycamore {Acer pseudoplatanus ) 


Forced flexural vibrations 


Tropical species (mean on 54 
species) 


Brancheriau et al. 

( 2010 ) 


Free vibration in bending 
and compression 


Table 2 Nominal dimensions and nominal properties of the tested timber beams 


Strength Fength Cross-section 
class 


Mean density Mean modulus of Mean shear 

(kg/m 3 ) 


(m 


(m x m) 


elasticity (GPa) modulus (GPa) 


Solid wood C24 
Glulam 


6 


0.220 x 0.070 420 
0.404 x 0.088 490 


11.0 


0.69 


GF32c 6 


13.7 


0.78 


Temperature and relative humidity were not controlled, but were assumed to 


correspond to common indoor conditions, respectively, 20 °C and 40 %. Each 
timber beam was simply supported with a symmetric overhang. Supports used were 
constructed of either rigid steel tripods or sections of thick steel cylinders. Teflon 
sheets were added in between the timber beam, and the steel supports in order to 
minimize friction and other sources of structural damping. The unique combination 
of span, mode, and orientation is referred to as a configuration. The configuration 
nomenclature is used as follows. The different test configurations are designated by: 

SW 


Solid Wood; GL 


Glulam; E = Edgewise orientation; F = Flatwise 
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orientation; second last number 
“SW-E-3-1” corresponds to the evaluation of the damping ratio £ for the first mode 
of transversal vibration for a solid wood beam lying edgewise with a span of 3 m. 


span; last number 


mode number. For instance, 


Experimental evaluation of damping ratios, fundamental frequencies, and mode 
shapes 


The modal hammer “heavy duty type 8208” from Briiel & Kjaer was used to set the 
beam into motion. A soft tip was employed in order to excite lower frequencies. 
Transient vibrations due to modal hammer impact were recorded by one ceramic/ 
quartz impedance head Kistler accelerometer type 8770A50 screwed into the beam. 
The load and acceleration time series were then digitalized and processed by a 
dynamic signal analyzer. The sampling frequency was fixed to 1,000 Hz, and 5 s data 
were recorded for each impact. An experimental modal analysis software was 
provided by National Instruments (2011) to record and process the data, using the 
graphical development environment Lab VIEW. Experimental modal analysis (Ewins 
2000) was used for determining the fundamental frequencies, the damping ratios, and 
the mode shapes of the timber beams, with the fundamental assumption of small 
damping. The frequency response function //, given with respect to the circular 
frequency co , uses the fast Fourier transform (FFT) to relate the input signal spectrum 
F(co) from the hammer and the output signal spectrum X(co) from the accelerometer as: 


X(co) 


H(co) 


2 


F(co) 


A linear average H 

impacts was performed for each evaluation. Identification of transfer function models 
was performed by curve fitting the averaged frequency response function. The 
parameter identification method was based on the Frequency Domain Direct 
Parameter Identification fitting method, which is a frequency domain multiple 
degree-of-freedom modal analysis method suitable for narrow frequency band and 
well-separated modes. The experimental setup is displayed in Fig. 1 and is similar to 
the one used to predict the presence of decay in logs (Ouis 2000). The present method 
is considered nondestructive (Ouis 1999) since the hammer impact is soft enough not 
to inflict any damage to the beam or modify its properties. This also allows an 
unlimited number of repeated measurements to be performed on each specimen. 

The driving point method (De Silva 2005) was used to evaluate the damping 
ratios and the fundamental frequencies. The impact and the data recording took 
place at the same location, 2.5 m from one end of the beams. All beams were 
evaluated for three different spans (3, 4, and 5 m) and on their four different faces 
A, B, C and D, where A and C represent the flatwise orientation, while B and D 
represent the edgewise orientation, as illustrated in Fig. 2. 

The roving hammer method (De Silva 2005) was used for mode shape 
measurements, in order to extract frequency response functions at several locations 
along the beam. In that case, 13 impacts were performed consecutively, every 0.5 m 
along the timber beam. The accelerometer remained screwed 2.5 m from one end of 


of the estimated frequency response function over three 


average 
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Fig. 1 Experimental setup 


e 


(b 


Fig. 2 Beam faces and orientations: a Face A and C, flatwise orientation, b Face B and D, edgewise 
orientation 


the beam during all consecutive impacts. The experimental evaluation of mode 
shapes was performed for one solid wood beam only. Three spans: 3, 4, and 5 m 
were considered. The first four modes were estimated for each span of the flatwise 
orientation, while only the first mode was recorded for each span of the edgewise 
orientation. 


General statistical methods 


Statistical indicators 


Various statistical indicators were evaluated for different sets of data X. In 
particular, the estimated mean value X, the standard deviation Sx , the median m x , 
and the coefficient of variation COV^ were evaluated. In addition, both the 95 % 
confidence intervals for the mean CI X and for the standard deviation CIs x were 
evaluated. They are, respectively, defined as: 


S x 


S x 


Cl 


X 


; X + 4).025 


^0.025 


X 


3 


(n~ 1) 


(n 


1) 


CIsv 


S x : 


S x 


2 


2 


%.025 


%.975 


where t om s 
the right, n 


t value with n 
number of evaluations, X 0025 an d X 2 


1 degrees of freedom, leaving an area of 0.025 to 

chi-square values with 


0.975 
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1 degrees of freedom, leaving areas of 0.025 and 0.975, respectively, to the 
right. The normalized length of the 95 % confidence interval for the mean X x is 
defined as: 


n 


2 


Sx 


(4) 


2 


^ 0.025 






For a normal distribution, the 99th percentile is defined by: 

P(X > X 99 ) = 1 % 

with X 99 = X + 233S X 

The Pearson product-moment correlation coefficient estimate r (Walpole et al. 
2007) was used to measure the strength of linear dependence between quantities. 
Similarly, the modal assurance criterion matrix (MAC) was used to measure the 
linear consistency (Allemang 2003) between two vectors u and e\ 


5 


2 


({uiV{ej}) 


( 6 ) 


MAC,, 


U 


{{utY {ui})({ej} T {ej}\ 


Statistical tests 


Statistical tests were also performed with a p value approach, associated with a level 
of risk of 5 %. The Anderson-Darling test was used to test the statistical hypothesis 
of normality, with the following hypotheses: 

H 9 : the data follow a normal distribution 
H i : the data do not follow a normal distribution 


7 


A one-way ANOVA test was performed in order to test the statistical hypothesis of 
equal means, which involves the following hypotheses: 

Ho : the means are equal 
Ho : at least two of the means are not equal 

The Leven/Brown-Forsythe test was performed in order to test the statistical 
hypothesis of equal standard deviations, and the following hypotheses were 
considered: 


( 8 ) 


Ho : the standard deviations are equal 

Ho : at least two of the standard deviations are not equal 


9 


Design of experiments approach 


Statistical Software (2010) was used to capture differences between 
response readings (outputs) for different groups of the input changes. For each 
factor, a number of levels are defined to represent the range of interest. Differences 
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<£) Springer 









Wood Sci Technol (2013) 47:1033-1050 


1040 


are then attributed to the factors acting alone (called a single effect) or in 
combination with another factor (called an interaction). 


Statistical treatment of damping ratio evaluations 


The probability distribution for all available evaluations for each configuration was 
assumed to follow a normal distribution. The process of detecting and treating data 
to be discarded is adapted from recommended practice published by the Ministry of 
Environment in British Columbia (CA) (Ministry of Environment of British 

Columbia 2001). 

Evaluations in the upper 1 % of the distribution were deemed unlikely to belong 
to the assumed population. The 99th percentile defined in Eq. (5) was therefore 
used as a lower boundary. An additional gap check was performed to investigate 
the inconsistency of the suspected evaluations. The associated statistical test 
compares the distance between the suspected evaluations and the next non- 
suspected data point to the predicted distance between the two largest values of the 
considered sample. These two largest values of the sample were reasonably 
assumed to correspond to the following percentiles of the assumed distribu¬ 
tion: 100 (A — 1 )/N and 100 (N — 2)/N. Any evaluation showing a distance to the 
next data point larger than twice the predicted distance was thus discarded as 
inconsistent. 


Repeatability and reproducibility investigations 


Repeatability is defined (International Organization for Standardization 2010) as: 
“the precision under observation conditions where independent test/measurement 
results are obtained with the same method on identical test/measurement items in 
the same test or measuring facility by the same operator using the same equipment 
within short interval of time.” Reproducibility is defined (International Organization 
for Standardization 2010) as: “the precision under observation conditions where 
independent test/measurement results are obtained with the same method on 
identical test/measurement items in different test or measurement facilities with 
different operators using different equipment.” 

In the present study, the repeatability and reproducibility were investigated by 
evaluating the damping ratio £ for the first mode of a single solid wood beam in the 
configuration SW-F-5-1, through a series of 10 evaluations, with 3 averaged hits per 
evaluation. The test series were performed by 10 different operators, 6 males and 4 
females with different weights and heights. Each operator was asked to perform 30 
hits using the same equipment, same test/measurement items, and same procedure. 
Only one operator, operator J, had previous knowledge of the equipment and the test 
procedure. 

Repeatability was investigated by evaluating the normality of the evaluations 
performed by each operator by means of the Anderson-Darling test defined in 
Eq. (7). Reproducibility was investigated by comparing the means and standard 
deviation of the series obtained for each operator, by means of a one-way ANOVA 
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Table 3 Selected factors and levels for the design of experiments 


Levels 


Factors _ 

Accelerometer location 


X2 


XI 


Impact location 


X2 


XI 


Span 


3 m 


4 m 


5 m 


Beam orientation 


Edgewise 


Flatwise 


hypothesis test defined in Eq. (8) and the Leven/Brown-Forsythe hypothesis test 
defined in Eq. (9), respectively. 


Parametric study 


A parametric study was performed on a single solid wood beam following a full 
factorial design of the experiment approach. The damping ratio £ for the first mode 
was evaluated for different combinations of spans, orientations, impact loading 
locations, and acceleration recording locations. Studied factors and their respective 
levels are reported in Table 3. The flatwise orientation results and the edgewise 
orientation results were computed as the mean values of evaluations related to faces 
A and C, and related to faces B and D, respectively. Locations XI and X2 
correspond to a distance of 2.5 and 3 m from one end-section of the beam, 
respectively. 


Numerical evaluation of mode shapes 


The commercial software Abaqus (Abaqus Analysis User’s manual, version 6.9 
2010) was used to perform the numerical analyses. Both solid wood and glulam 
beams were modeled using 4,536 quadratic continuum elements with reduced 
integration. The mesh was chosen with at least 9 elements to compose the height of 
the beam, and at least 8 elements along the width of the beam. Approximately, 60 
elements were used along the beam in the longitudinal direction. The nominal 
material properties given in Table 2 were used. 


Results and discussion 


Validation of the experimental protocol 


Investigations on repeatability and reproducibility were performed using a panel of 
operators, as described in “Repeatability and reproducibility investigations,” in 
order to evaluate the consistency and robustness of the impact testing method. Out 
of 100 records, one acquisition could not be processed for technical reasons. The 
results are reported in Fig. 3 and emphasize the presence of an outlier for operator 


E. 


The p value for the Anderson-Darling test was computed as 0.310 for a level of 
risk of 5 %, which reveals that the normal distribution has a significant fit to the 
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Fig. 3 Box plot describing all damping ratio evaluations for all operators 


Damping ratio 


Fig. 4 Damping ratios histogram and normal distribution 


data, as shown in Fig. 4. Except operator E, whose higher coefficient of variation 
(11 %) is due to an outlier, all operators exhibit a low coefficient of variation similar 
to that of operator J (3 %), which demonstrates that experience is not needed to 
produce good repeatability. 

According to the one-way ANOVA test, the mean values show no significant 
differences (p value = 0.162) from one operator to another one. According to the 
Leven/Brown-Forsythe test, the standard deviations (p value = 0.062) have no 
significant difference from one operator to another. Since the experimental protocol 
demonstrates good reproducibility, the damping ratio evaluations are therefore 
independent of the strength and skills of the operator. 

Excellent agreement between damped experimental fundamental frequencies and 
undamped numerical fundamental frequencies is observed in Fig. 5. In addition, 
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Fig. 5 Comparison of experimental and numerical natural frequencies 


large diagonal values (~0.98) in the MAC matrix defined in Eq. (6) reveal good 
agreement between damped experimental mode shapes and undamped numerical 
mode shapes for each fundamental frequency. The assumption of small damping, 
which is mandatory for applying experimental modal analysis theory, is therefore 
valid. 


Eventually, neither impact locations (modal hammer) nor the measurement 
locations (accelerometer) were found to have a significant influence on the damping 
ratio £ in the parametric study described in “Parametric study” section. This 
indicates reciprocity, which is a fundamental characteristic of experimental modal 
analysis methods. 


Reliability of the experimental evaluations 


The statistical treatment of damping ratio evaluations described in “Statistical 
treatment of damping ratio evaluations” section was performed twice. In total, 14 
evaluations were discarded out of 420 evaluations: three evaluations for glulam 
beams and 11 evaluations for solid wood beams. The low percentage of discarded 
values: 


3 % indicates consistent results. In addition, the obtained experimental 
damping evaluations compare well with the results from previous studies reported in 
Table 1, although slightly higher. Most of the measurements reported in Table 1 
were obtained from small-scale specimens with assumed low inhomogeneity, while 
the present study used specimens of structural dimensions, thus exhibiting larger 
statistical inhomogeneity. The larger statistical inhomogeneity is most probably one 
of the main reasons for the slight increase of damping. 

The mean values of damping ratio evaluations \ and of natural frequencies / per 
configuration are presented in Table 4, together with number of evaluations, number 
of discarded evaluations, and other statistical indicators. In general, the number of 


42 Springer 









Wood Sci Technol (2013) 47:1033-1050 


1044 




o 


c/5 OD G 

HAS 


coin^j-coocvocoocin 

(N(Nr"-VOVOOCcOOC(N 

ov i> © co co in 


voinoc^^tovr-co 


On co 
co vo 


0-) • 


03 


73 


co 


in 


fl c3 > 

< Q a. 


On 


CO 


CO 


in 


in 




ovint^vovoovocvoin 

CO (N 


OC C\ CO (N CO 

O (N 7f CO (N 


(N CO (N CO r- 
CO <N 


(N 


(N 


r- 


KJ> 




> 


vo ri ci cc 
rn h oo 


-xt vn vn ri vn a co 

cs vo cn tj-" vd vn i> tj-" 

(N (N 


CN 




oo o 


o 




ON CO OO Tt 
(N <N 


OO 


co co 

<N <N co co 


u 


(N 




(N 




in oo Tt oo 
in oc 


tj- in co oc 

co vo ov co 


oc o\ 

CO ^ C h 


h' ^ OC If) ^ 


(N 


(N 




OC GV (N (N CO 

oc (N in 


(N t 00 CO 

(N (N 


in vo 


(N ^ in oo 
(N in ^t 


r- 


CO CO 
VO CO (N 


7i- 


CO oc 


(N 


rs 


rs 


#N 


rs 


#N 


rs 




rs 


rs 


#N 


rs 


#N 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


OO h VO CO (N 


t"^ <N "\f GV 


(N (N OV 

(N (N 


h oc vo ^ 


CO oc 


in 


CO 


CO 




CO 


o 


VO <N ON ON (N in 
04 04 


04 ON co ON co 


On OO 


04 04 CO 


co 


co 04 


04 


04 


OO 


VO) 


CO 


Hjl 


to 


^j-ovovininvoovocco^t 

CO CO CO oc oc 


co oc r- vo oc vo 

in in 


GV OC 

co in r- 


vo 


oc 


(N 


(N 


(N 


CO 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


rs 


VO 


VO CO (N 

ov co in oc co 


vo in vo 
vo vo 


hochco^roh^ 
co in h h 


GV 


OC CO 


r- r- 


l*U» 


u 


in ^j- ^t 
co vo 


r^^j-invoocco^|-(NOvco(NinGNco 

vocoinvoGvocvo^tov^tinin 


in in gv 
in gv r- 


CO 


(N 


(N 


(N 


(N 




N 


X 


fl 


tvovoott 
vo 1- 1- co in 


(N in in (N ov co 
CO ^ vo 


OC h G) VO 

CO Tt (N 


vo 


GV 


O 


(N 


VO 


co r- 


ctf 




p 




VR 


cn 




fi 


fi 


O 


73 


O 


O 


G . 




o 


G 73 


cd 




P P 


G 


a 


a 


o 


cd 


3 C/3 

2; ^ 


> 


C/5 


(N 


G 


O 


cd 




C/3 


6D c^H 


o P 


fl 


O 


Oh 


G 


a 


cd 


p 


cd 


a 


73 


cd 


P 


OC xj- G>. oc 


(N (N 


m(NOCcOCOGVGV(Nr^ 


vo in 


> 


73 


Z 


(N 


(N 


(N (N 


(N 


G 


G 


cd 


P 


cd 


C 


> o 


W 


(N 


(N 


(N 


cd 


(N 


(N 


(N CO 


(N CO Tt 


(N CO 


i 


i 


i 


i 


i 




i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i a 


CO CO VO VO CO CO 




I 


I 


I 


I 


I 


6D 


i 


i 


i 


i 


i 


i 


i 


I I I I I I I I I I I I I I I I hJ hJ hJ kJ kJ kJ kJ 

oooooooooooooooo£££££££ 




y=! 


pfi 


c 




o 


H U 


^ Springer 







Wood Sci Technol (2013) 47:1033-1050 


1045 




o 


H fl 3 


00 oo 


in 


<30 • 


03 


m 


o\ 






c c3 > 

< Q ^ 


m <n ^r 


m 


m 




Tf On vo 
(N VC OC 


co ^|- 
(N (N 


OO 




co 




<N 




> 


oo r- (N 
vn in 


in o cn in 

o\ t> 


O 




in 


U 


(N 


vo in 

r-inCNCN'sOi^inr- 




oo 






(N 0> 


OO OO h O 

C\ (N (N 


G> 


co 


(N co 


co 


• Cn 


rs 


#N 


rs 


#N 


rs 


rs 


rs 


rs 


t m m m t 
<N in m 




oo 






CO 


U 


oo in in (N 


in in On co in co 

CO 00 vo 


00 ^1" 


(N 




to 


^t 


in 


CO 


00 oo VO G> 


G> 


oo 


On 


(N 


CO 


• cn 


• CN 


• #N 


• CN 


• #N 


• CS 


• CN 


• CN 


• CN 


oo r- 

VO ^t 


^ ^ in r- 

vo o m vo 


r- 


r- 


|H_P 


u 


co^t^sor-cNcNooaN 

in r- r- in r- 


CO 


oo 


(N 


(N 




N 


X 


vo t> t> 
(N co in 


in o\ 
(N in oo 


CO 00 


C/3 




c 


o 


73 


o 


oo . 




CD 73 




P 3 


B 


O 


03 


3 


> 


(N 


CO 


<30 


<+H 22 

O C 


O 


<30 


o3 


3 


a 


P 


P 


P" (N 


On VO CO 


in vo in 


> 




(N 


(N 


73 


<30 


oo 


p 


C 


C 


fl 


o 


o 


o 


CO Tf 

coco^tp-p-p-inininin 

. 

H,tL,tL,tL,tL,tL,tL,tL,tL,tL, 


(N cn ^r 


(N m ^t 


cd 


5-1 


I 


I 


I 


I 


I 


I 


I 


I 


7f 


p 






<-C 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 


pfi 


fl 






O 


H U 


GOGOGOGOGOGOOOGOOOGO 


ringer 






Wood Sci Technol (2013) 47:1033-1050 


1046 


damping evaluations decreases when the span is decreased and/or when the mode 
number is increased. Higher modes were difficult to obtain for shorter spans because 
the modal hammer supplied only a limited amount of energy to the structure. As 
there were few available evaluations, this led to wider 95 % confidence intervals 
and to a larger 2«*. Three configurations out of the 33 are based on three or less 
evaluations: GL-S-3-1, GL-W-4-4, and SW-W-3-4. In addition, four configurations 
exhibit a 2^ larger than 100 %: GL-S-5-2, SW-S-3-2, SW-W-3-3, and SW-W-4-4. 
Results related to these specific configurations may therefore not be reliable enough 
for further use. 

However, results concerning long spans, low modes, or flatwise orientations 
revealed narrow confidence intervals and low values of S^ and of COV which 
indicates good reliability. Similarly, the observed low values of 2^ indicate a limited 
scatter of the evaluated damping ratios, a large number of evaluations, or both. In 
addition, large p values for the Anderson-Darling test indicate consistent 
evaluations for most configurations. General good reliability of the obtained 
experimental damping evaluations is therefore demonstrated. 


General trends 


Difference glulam/solid wood 


Although the total amount of damping evaluations is similar for glulam beams and 
solid wood beams, glulam beams show generally smaller COV^ and smaller 2^ than 
solid wood beams for most of the configurations. This is in accordance with the 
measured densities, reported in “Experimental setup” section. This confirms a 
generally observed lower scatter in the natural properties of glulam compared to 
solid wood. 


Shear deformation effect 


The sensitivity of damping ratio evaluations to experimental conditions was 
analyzed using the systematic design of experiments approach, as defined in 
“Design of experiments approach” section. The governing effects were found to be: 
(1) the beam orientation, (2) the span, and (3) the interaction between the beam 
orientation and the span. From the results reported in Table 4, the evaluated 
damping ratio appears generally larger for edgewise orientation, for higher modes, 
and for shorter spans. In particular, the damping ratio seems greatly dependent on 
the mode number, and this dependence increases as the span decreases. Glulam 
beams, with a cross-section that has a height-to-width ratio larger than that of solid 
wood beams, exhibit higher evaluated damping ratios (+20 %) than those of solid 
wood beams. The reported influences are captured by the influence of shear 
deformation, which is hence observed to significantly increase damping. Similar 
findings were reported concerning Timoshenko wooden beams (Nakao et al. 1985) 
and laminated composites (Hwang and Gibson 1991). 
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Experimental investigation, limitations, and uncertainties 


In the present study, uncertainties are related to the measurement technique, the 
processing technique, the material itself, and the specimen geometry. The known 
limitations of the experimental setup were initially accounted for. In order to ensure 
that the evaluated damping was representative of only material damping, Teflon 
sheets were used to reduce friction and hence the structural damping contribution. 
Supports were resting on isolated massive concrete blocks, so that the least amount 
of vibration was transmitted to the supporting structure or to the surrounding fluid, 
in order to minimize the contribution of fluid damping. 

In addition, great care was taken when sampling and windowing the signal in 
order to reduce signal distortion for the estimated frequency response functions. 
When processing the recorded data to extract modal properties, curve fitting was 
used for parameter identification and the residual errors were assumed to remain 
small. 

Scatter may be expected for all wood material properties because wood exhibits 
high variability. No quantification of variability of material damping is available in 
standards, but the natural variability of the longitudinal elastic modulus of GL32c 
(European Committee for Standardization 1999) and of C24 (European Committee 
for Standardization 2009) is quantified by a coefficient of variation of 12 and 20 %, 
respectively. Variability of the material may also induce uncertainties in the 
specimen geometry. Here, differences were observed for pairs of opposite faces and 
were probably due to initial warping. Larger differences were observed in solid 
wood beams, whose shapes are less stable than those of glulam beams. Foster 
(1992) made similar observations and concluded on the variations of the beam 
surface from one face to the opposite one. 


Correlation investigations 


The correlation between damping, natural frequency, and density was evaluated by 
means of the Pearson product-moment correlation coefficient, as defined in 
“Statistical indicators” section. In general, medium correlation (r = 0.46 for Solid 
Wood, r = 0.15 for Glulam) was observed between density and fundamental 
frequency. The longitudinal modulus of elasticity is known to be medium correlated 
with the density (Joint Committee on Structural Safety 2007), so it is expected that 
an increase of density induces an increase of stiffness and thus an increase of the 
natural frequency. Poor correlation (r 
Glulam) between damping ratio and density demonstrates that the evaluated 
damping ratio is not significantly influenced by the density. Similar findings were 
reported by Obataya et al. (2000), who found a high negative correlation between 
damping and specific modulus—that is, ratio modulus of elasticity/density. They 
finally assessed that both quantities depended on the microfibril angle but little on 
density. Poor correlation (r 
also observed between damping ratio and fundamental frequency. This is in 
accordance with previous findings (Krueger and Rohloff 1938; Ono and Norimoto 
1985; Nakao et al. 1985; Foster 1992). Even if wood is considered a viscoelastic 


0.01 for Solid Wood, r 


0.11 for 


0.04 for Solid Wood, r 


0.02 for Glulam) was 
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material (Ouis 2002) whose properties depend on frequency, the frequency range 
per configuration is very narrow (less than 5 Hz). The frequency dependence is 
therefore unlikely to be observed in general. 

However, when computing the Pearson product-moment correlation coefficient 
for a group of configurations, rather than a single configuration, correlation can be as 
high as excellent between damping ratio and fundamental frequency. For instance, 
r > 0.9 for glulam beams in flatwise orientation, for any span. This indicates strong 
mode number dependence of the material damping for a given system, that is, when 
material, orientation and span are fixed. Previous studies (Lazan 1968; Kume et al. 
1982) also underlined the specific influence of the mode shapes. 


Influence of shear deformation and mode shape 


Experimental mode shapes are observed to be similar for configurations presenting 
identical span and mode number but different orientations, as long as the shear 
deformation is not significant compared to bending deformation (i.e., for the first 
modes and for longer spans). As expected, similarities decrease for higher modes 
and shorter spans. 

For slender Euler-Bernoulli beams, whose plane sections are expected to remain 
plane after deformation, the shear V EB is related to the flexural rigidity EL For 
Timoshenko beams with a smaller length-to-height ratio, the deformation due to 
shear is considered. Both types of shear are expressed with respect to the mode 
shape: y = f(x ): 


0/C 


/// 


El 


Ely 


Feb 


( 10 ) 


dx 


kAGy' 


V T 


curvature, k = shape factor, A 


cross-section area and G 


shear 


where k 
modulus. 

Consequently, the shear V depends on specific characteristics of the mode shape: 
its third derivative y'" and/or its slope /. Since it is difficult to quantify the shear 
directly from the experimental measurements, it is therefore possible and more 
convenient to capture the material damping variations via the influence of the mode 
shape and its derivatives. 


Conclusion 


The present study focused on examining the material damping in timber beams of 
structural dimensions for different experimental configurations. A total of 22 timber 
beams of two different types were subjected to flexural vibrations through the 
impact test method. Damping evaluations were performed for various configura¬ 
tions, which included different spans as well as orientations (edgewise and flatwise). 
A total of 420 evaluations were performed, out of which 14 were discarded 
following a rigorous statistical process. Statistical indicators to the damping 
evaluations were provided together with their mean values for each configuration. 
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The initial assumption of small damping was validated by means of mode shape 
analysis. The reliability of the results was assessed by carefully evaluating the 
consistency of the data and by concluding that no significant differences were due to 
the operator and/or his/her skills. In addition, the natural variability of wood and the 
limitations of the experimental setup were thoroughly discussed in qualitative and 
quantitative manner. 

General trends were that the evaluated damping ratio increases with higher 
modes, shorter spans, and the edgewise orientation as compared to the flatwise 
orientation. The investigated influence of density on the evaluated damping ratio 
found no significant correlation. Correlation analysis between the evaluated 
damping ratio and the corresponding natural frequency revealed that the mode 
number is a significant parameter for a given system, but not for a given 
configuration. Instead, the shear deformation was found to be the governing factor 
to explain the variation of the damping ratio from one configuration to another. 
Shear deformation was finally conveniently evaluated by mode shape 
characteristics. 

The comprehensive material damping database presented in this study can be 
used in further investigations to enhance general knowledge on damping mecha¬ 
nisms and improve the prediction of comfort properties for timber structures in 
particular. 
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